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1. Introduction 



In this paper we study the low energy dynamics of supersymmetric QCD (SQCD) in 
the presence of certain F-term deformations. The starting point of our analysis is an N = 1 
supersymmetric gauge theory with gauge group SU (N c ) and Nf > N c flavors of chiral 
superfields in the fundamental representation of the gauge group, Qf , Q l a (a = 1, • • • , N c ; 
i = 1, • • • , Nf). This theory has a global symmetry 

SU(N f ) x SU(N f ) x U(1) B x U(1) R (1.1) 

and a non-trivial moduli space of vacua, which has been extensively studied and is rather 
well understood; see e.g. ]l|-|6| for reviews. 

A natural question is what happens when we deform the theory by adding a general 
superpotential that preserves a particular subgroup of the global symmetry (|1.1|), such as 
the non-chiral subgroup SU (Nf)^ ag xU(1)b- A class of superpotentials with this property 
is 

W el = Y -m n TtM n (1.2) 
^— ' n! 

n=l 

where the meson field 

M) = Q l Q 3 (1.3) 

is an Nf x Nf matrix (the color indices are summed over in ( |1.3| )). Terms with n > 1 
in ( |1 . 2|) are non-renormalizable, which is reflected in the fact that the couplings m n have 
dimension [m n ] = 3 — 2n at the free fixed point. One can think of the superpotential ( |1 . 2j ) 
as providing an effective description below a certain energy scale. 

In the case m n = mi5 n ,i ( |1.2|) is a mass term for Q, Q. It has been known for a 
long time that the resulting theory has N c supersymmetric vacua, in accordance with the 
Witten index. More recently, it was found that for N c < Nf < |iV c and small mi there 
are metastable non-supersymmetric ground states as well. Such states might be useful for 
describing supersymmetry breaking in nature. 

The purpose of this paper is to study more general superpotentials of the form ( |1.2| ) . 
We will mainly discuss the case where only the two lowest terms in ( |1.2| ) are non-zero, i.e. 

I 

W el = miTrM + -m 2 TrM 2 (1.4) 

i 



but will also comment on higher order superpotentials. We will see that such superpo- 
tentials lead generically to a rich landscape of supersymmetric and non-supersymmetric 
vacua, and explore some of their properties. Related works include 0-0. 

Like in 0, we will find it useful to utilize the Seiberg dual description of SQCD, in 
which the gauge group is SU(Nf — N c ), and the meson (|1.3| ) becomes a gauge singlet 
field. We will analyze the dynamics in both the electric and the magnetic descriptions and 
compare them. 

The deformed SQCD with superpotential ( |1.4| ) has a simple embedding in string the- 
ory, along the lines of |TI||50[]. In a companion paper |JT| we describe the relevant string 
construction, and in particular its connection to the gauge theory results of this paper. 
We find that the brane description provides a complementary picture to the gauge theory 
one. 

The plan of the paper is as follows. In section 2 we construct the supersymmetric 
ground states of the theory ( P-.4Q , and verify that the electric and magnetic descriptions 
give rise to the same vacuum structure once all the relevant quantum effects have been 
included. In section 3 we describe metastable states in this theory. In section 4 we discuss 
our results and comment on generalizations. 

2. Vacuum structure of deformed SQCD 

2.1. A first look at the vacuum structure 

The N = 1 supersymmetric Yang-Mills theory with gauge group SU(N C ) discussed in 
the previous section, which we will refer to as the electric theory, is equivalent in the infrared 
||32j| to another gauge theory, which we will refer to as the magnetic theory. The latter has 
gauge group SU (Nf — N c ) and the following set of chiral superfields: Nf fundamentals of 
the gauge group, q % , q~i, i = 1, ■ • • , Nf, and a gauge singlet Mj which, as suggested by the 
notation, is identified with the gauge invariant meson field (|1.3| ) in the electric theory. 

The magnetic quarks q, q and meson M are coupled via the superpotential 

W mag = jqiAQq* . (2.1) 

The scale A is related to the dynamically generated scales of the electric and magnetic 
theories, A e , A m , by the scale matching relation 0] 

^ c -N fA HN f -N c )-N f = { _ )Nf -N cA N f {22) 
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We are interested in adding to the electric Lagrangian the superpotential 



W ei = -Tr(QQ) 2 - mTr(QQ) = — TrM 2 - mTrM (2.3) 
2 2 



which has the form (|1.4j ), with mi = — m and mi = a. In the magnetic description this 
corresponds to deforming ( |2.1j ) to 

WW = \qM)q J + 7^ TrM 2 - mTrM . (2.4) 



Since the superpotential ( |2.4| ) is quadratic in M, we can integrate this field out. The 
resulting superpotential for the magnetic quarks is given by: 



^Tr(qq) 2 - mTr(qq) 



(2.5) 



Comparing to (|2.3|) we see that the magnetic superpotential has the same qualitative form 
as the electric one. 

Conversely, one can write the electric superpotential ( |2.3| ) in a way similar to the 
magnetic one by integrating in a gauge singlet field N: 

W e i = -^-Q*N?Qj - y TriV 2 + m e TrN . (2.6) 

Requiring that integrating out iV leads back to (|2.3| ) gives 

1 m e 
a = — , m = — - . (2.7) 

a e A 2 a e A 

Comparing ( p.6| ) to (|2.4j) we see that the two superpotentials have very similar forms. Note 
also that we chose the normalization of the field iV (which transforms as a singlet under 
the electric gauge group) such that it is identified in the infrared with the magnetic quark 
bilinear N = qq. This is the dual version of the relation between the singlet meson M in 
the magnetic theory and the bilinear in electric quarks, ( |1.3j ). 

Although the two forms of the magnetic superpotential, ( |2.4|) and (|2.5|), describe the 
same long distance physics, there is a physical difference between them, which becomes 
important for small a. In this limit, the mass of M in ( |2.4j ) is small, and at energies above 
that mass this field should be included in the dynamics. The description ( p.5|) does not 
contain this field; it coincides with ( |2.4j ) only at energies well below the mass of M. Similar 
comments apply to the electric superpotentials (|2.3|), ( [2.6| ) in the limit a e — > 0. 



To recapitulate, we see that the analysis of the vacuum structure in the electric and 
magnetic descriptions is essentially identical, up to the replacement iV c <-> Nf—N c , a <-> a e , 
m «-> m e , etc. We will verify later that the two descriptions lead to the same vacuum 
structure, in agreement with Seiberg duality. 

At first sight, it actually seems that the electric and magnetic theories have different 
vacuum structures. Consider, for example, the classical magnetic superpotential (|2.4j). To 
find the supersymmetric vacua we need to solve the F-term constraints 

My =o , 

QiMj =0 , (2.8) 
— qiq-' =md J i — olM{ . 

From (|2.8j ) we learn that M satisfies the matrix equation 

mM = aM 2 . (2.9) 



On solutions of the equations of motion one can choose M to be diagonal. Equation 
implies that its eigenvalues can take only two values, and — . Without loss of generality 
we can take 

M=(° ° ) (2.10) 



where k = 0, 1, 2, • • • , Nf and I n is the n x n identity matrix. Plugging (|2.10|) into the last 
line of ( |2.8| ) implies that 

( mA Q Ik I) • (2.11) 



qq 



The rank of the matrix on the left hand side is at most Nf — N c . Hence, one must have 
k<N f - N c . 

Note that the vacua ( |2.10|) , ( |2.11| ) go off to infinity in field space as the deformation 
parameter a (|2.4j ) goes to zero. In particular, for a = supersymmetry is spontaneously 
broken. This is consistent with the fact that for a = the classical superpotential ( |2.4| ) 
has an unbroken U(1)r symmetry, with R q = Rq = 0, Rm = 2. For a^0 this symmetry 
is explicitly broken, and one expects [[J3| that supersymmetric vacua exist. 

Expanding around the solution (|2.10| ) , ( |2.11|) one finds that the only massless degrees 



of freedom are gauge fields and fermions associated with pure N = 1 SYM corresponding to 
the unbroken gauge group SU (Nf — N c — k). As mentioned above, quantum mechanically 
this theory has Nf — N c — k vacua with a mass gap (for Nf — N c — k > 2). Thus, the 



above analysis implies that the theory with superpotential ( |2.4|) has the following number 
of vacua (up to Nambu-Goldstone bosons associated with broken global symmetries): 



N f -N c -1 



N mag 



1+ J2 (N f -N c -k) = l + hN f -N c )(N f -N c + l) . (2.12) 



s=0 



The 1 in (|2.12| ) corresponds to the case k = Nf — N c where the expectation values of 
the baryon fields b = q N f~ Nc and b = q N f~ Nc are non-zero, which we will refer to as the 
baryonic branch. 

As mentioned above, the analog of the magnetic superpotential (|2.4| ) in the electric 
theory is (|2.6|). This superpotential clearly leads to the same vacuum structure as ( |2.12| ) 
with the replacement Nf — N c — > N c , 



N r -1 

k=0 



Ne-l -, 

1 + ( N ° - k ) = 1 + oN c (N c + 1) . (2.13) 



Thus, the electric and magnetic answers are in general different, whereas Seiberg duality 
implies that they must agree. The resolution has to do with quantum corrections to the 
superpotentials of the electric and magnetic theories. We next turn to a discussion of these 
corrections, first in the special case Nf = N c + 1, and then in general. We will see that 
after including quantum effects the number of supersymmetric vacua is given by 

iV vac = max(iV el , iV mag ) (2.14) 

in both the electric and the magnetic theories. 

2.2. Quantum corrections for Nf = N c + 1 

This case is particularly simple, since the magnetic gauge group is empty. At the same 
time, according to (|2.14j ) the number of vacua is in this case supposed to be given by the 
electric result ( p. 13 ) , which is larger than the magnetic one Q2.12 ). The magnetic quarks 



q 1 are proportional to the electric baryons B l , and it is convenient to use the latter as the 
fundamental degrees of freedom. The magnetic superpotential (|2.4D is known to receive 



an important correction proportional to detM, and takes the form |34 

1 
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lT liK1 , = ■ i2 ^ / _3 (BiM}B j - detM) + |TrM 2 - mTrM . (2.15) 



The new term in the superpotential leads to a correction to the F-term condition on the 
third line of ( |2.8| ), which now becomes 

BiB j - {detM)(M- 1 y i + A 2 e Nf ~ 3 (aM? - mSfj = . (2.16) 

To analyze the solutions of these equations we need to distinguish between the cases where 
the meson matrix M is regular and singular. Consider first the case where it is regular, 
which we will refer to as the mesonic branch. Then, the first two lines of ( |2.8|) (with 
q 1 oc B l ) imply that 

B i = B J = (2.17) 

while ( |2.16| ) takes the form 

Ae^'" 3 (aM 2 - mM) = (detM)I Nf . (2.18) 

As before, on solutions of the equations of motion we can diagonalize M. Since the left 
hand side of (|2.18|) is proportional to the identity matrix, the Nf eigenvalues of M must 



take at most two distinct values, which we will denote by x and y, and are related as 
follows: 

m 

x + y = — . 2.19 

a 

The vacua split into classes labeled by an integer < k < ^f, with k of the eigenvalues 
of M equal to y, and Nf — k equal to x. The upper bound on k takes into account the 
freedom of exchanging x and y. 

For k = 0, M is proportional to the identity matrix, M = xlN f - Plugging this form 
into (tnS|) we find 

Ae* /-3 (aa; - m) = x N f~ x . (2.20) 

There are Nf — 1 = N C solutions for x, and thus A^ c vacua with this form of M. For generic 
a, and in particular when a is sufficiently small, the A^ c solutions of ( |2.20| ) are all distinct 
and non-vanishing, as expected. 

For k 7^ the two distinct eigenvalues of M satisfy 

A 2 e Nf -\ax-m) =x N f- k ~ 1 y k , 

(2.21) 

Al N '- 3 (ay-m)=x N '- k y k - 1 . 



Using (|2.19|) we can rewrite both lines of (p.21| ) as 

-aAl Nf - 3 =x N '- k - 1 y k - 1 ■ (2-22) 
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Since y is linearly related to x, (|2.22|) is a polynomial equation for x of order Nf — 2. It 
has Nf — 2 = N c — 1 distinct solutions with x, y ^ (again, for generic a). Thus, for odd 
Nf the number of supersymmetric vacua of the magnetic theory is given by: 

iV c + \{Nf- 1)(N C - 1) = ^N C (N C + 1) . (2.23) 

For even Nf we have to analyze separately the case k = Nf / 2 for which the degeneracies 
of the two eigenvalues x and y are equal. In this case, due to the symmetry of interchange 
of x and y we actually have only \(Nf — 2) = \(N C — 1) distinct solutions of ( |2.22| ). The 
total number of vacua is in this case 

N c +\{N f - 2)(N C - 1) + \{N f - 2) = \N C (N C + 1) . (2.24) 

We see that for both even and odd Nf, the total number of vacua in the magnetic descrip- 
tion (|2.4|) agrees with ( [2.14|) , except for the contribution of the baryonic branch to which 
we turn next. 

So far we assumed that the meson matrix M is non-degenerate. The analysis needs 
to be modified when some of its eigenvalues vanish. It turns out that the only non-trivial 
case is the one in which exactly one eigenvalue of M vanishes. Thus, we take 

M = diag(M 1 ,M 2 ,---,M Nf ) (2.25) 

and consider the limit Mi — > with the rest of the Mj remaining finite. The first two lines 

of ( |2.8|) imply that the only non-zero components of 5, B are B , B\. Eq. ( |2.16| ) leads to: 

m 

Mj=-, j = 2,---,N f , 

a (2 26) 

B iB 1 = rr)*''- 1 +m Ar'- 3 . 



a 

Including this extra vacuum, which corresponds to the baryonic branch in the classical 
analysis, brings the total number of vacua into agreement with ( |2.14| ). 

If k > 1 eigenvalues of M go to zero, ( |2.16| ) has no solutions. Indeed, in that case the 
first two lines of ( |2.8|) imply that the non-zero components of B, B lie in the degenerate, 
k dimensional, subspace. In that subspace, (|2.16|) implies that 

oc 8) , (2.27) 

which is impossible to satisfy since the left hand side has rank 1 while the right hand side 
has rank k > 1. 

To summarize, for Nf = N c + 1 we conclude that the number of vacua of the magnetic 
theory is given by ( |2.14| ). The term that goes like detM in the magnetic superpotential 
( |2.15| ) is crucial for the analysis. The analog of this term for generic Nf > N c + 1 comes 
from quantum effects in the magnetic theory. We next turn to a discussion of these effects. 
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2. 3. Quantum corrections for Nf > N c + 1 

We would like to extend the analysis of the previous subsection to general Nf. Con- 
sider a vacuum in which the expectation value of the meson field M has r vanishing 
eigenvalues and Nf — r non- vanishing ones. As mentioned in subsection 2.1, the F-term 
constraint on the third line of (|2.8|) implies that r < Nf — N c . Therefore, at least N c 
eigenvalues of M must be non-zero. In particular, we can bring it to the form 

M=(M °) (2.28) 



M| 



where Mq is a non-degenerate (N c — 1) x (iV c — 1) matrix, and M is an (Nf — N c + 1) x 
(Nf — N c + 1) one, which may or may not be degenerate (but whose rank is at least one). 

Since M is non-degenerate, the flavors of quarks q, q that couple to it via ( |2.1|) are 
massive and can be integrated out at low energies. This leads to an SU (Nf — N c ) gauge 
theory with Nf — N c + 1 flavors q a , q& (a, a = 1, • • • , Nf — N c + 1), whose scale is given by 

A 2(N f -N c )-l = detMo A 3(Nf-N c )-N f ^ ^ 

This theory is of the sort discussed in the previous subsection. It can be described in terms 
of the gauge invariant observables 

N~ =q- a -q a , b a = q N f~ N ^ , b & = q N f~ N - . (2.30) 

The full magnetic superpotential has the form 

W mag = Tr (jMN + |M 2 - mM^j +Tr (|m o 2 - mM )+ (b ■ N ■ b — detiV 

_ L (2 ' 31) 

The field M appears quadratically in ( |2 . 3 1| ) and thus can be integrated out. The equation 

of motion of M sets it to 

M = —^(N - mAI Nf - Ne+1 ) . (2.32) 

Plugging this into ( |2.31| ) (and dropping a constant contribution to the superpotential) 
leads to 

Wuiag = "iv^ Qa n2 ~ mN ) + Tr (f M ° 2 - mM ° 



A 



(2-33) 



.,. x v . v (b ■ N ■ b - detN) 







The part of the superpotential that depends on N, b, b is very similar to the one analyzed 
in the previous subsection, ( |2.15|) . The new element is the dependence on Mq which needs 
to be taken into account. 

As before, the F-term equations of motion of N, 6, b have two kinds of solutions: one 
in which A" is non-degenerate, and another in which it has exactly one vanishing eigenvalue. 

Consider first the mesonic branch, where A" is non-degenerate. In this case it is 
convenient to go back to ( |2.31| ) and integrate out N, b, b. This amounts to setting 6 = 6 = 
in Q2.31|) and replacing A" by the solution of its equation of motion, 



-~- A^ detA" 

= A 3 W .- W .,-N, detM > -A-.+i ■ (2-34) 



Solving for N and substituting in ( |2.31| ) (using (|2.2| )) we get 



W mag = Tr (|M 2 - mM) - (N f - N c ) (jS^) (2-35) 
where M is the full meson matrix (|2.28|), which is non-degenerate in this case. Of course, 



the determinant term in ( |2.35| ) is nothing but the well known non-perturbative superpo- 
tential flfll for the meson field in SQCD. We could have gotten it directly from (|2.4D by 
assuming that M is a non-degenerate matrix and integrating out the massive quarks g, q. 
The fact that we got it with the correct coefficient is a check on the algebra. 

We can now generalize the discussion of the previous subsection and look for (non- 
singular) solutions of the F-term equations corresponding to the superpotential (|2.35| ): 



(detiW \ N f~ N c 
A 3N C -N f ) T Nf ■ (2-36) 

Again, the matrix M has only two distinct eigenvalues x, y satisfying the relation ( |2.19| ), 
and vacua are labeled by an integer k which keeps track of the number of times the 
eigenvalue y (say) appears. 

Considerations similar to those that led to ( |2.22|) give 



x 



N c -k y k-N f +N c = { _ a)Nf - NcA 3N c -N f (2 _ 37) 



Together with (gjg| ) this gives a polynomial equation for x, whose order depends on Nf, N c 
and k. 
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Consider first the case Nf < 2N C . For k < Nf — N c , (|2.37 ) is a polynomial of degree 
N c — k. Thus it has N c — k solutions. For Nf — N c < k < the degree of the polynomial 
and the number of solutions are given by 2N C — Nf. Summing over k one finds that the 
number of vacua (up to global symmetries) is ^N C (N C + 1), as before ( |2.23j ). 

For Nf > 2N C the picture is slightly different. For < k < N c one finds Nf ~ N c — k 
solutions, and for N c < k < Nf — 2N C solutions. The total in this case is |(iV/ — 
N c )(Nf — N c + 1), again in agreement with ( |2.14| ). 

It is not surprising that for Nf > 2N C the number of vacua agrees with the magnetic 
answer, while for Nf < 2N C it does not. The non-perturbative superpotential in ( |2.35| ) 



gives a correction to the classical superpotential for M, (|2.4| ), that goes like M N f~ Nc . For 
Nf > 2N C this correction is subleading at large M relative to the leading, M 2 , term, and 
one does not expect it to change the number of vacua. On the other hand, for Nf < 2N C 
it changes the behavior of the potential at infinity, and it is natural that the number of 
vacua changes. 

So far we assumed that the matrix N ( |2.30| ) is regular. As we saw in the previous 
subsection, the only other case we need to consider is the baryonic branch, in which the 
rank of N is Nf — N c (i.e. one of the eigenvalues of iV goes to zero). We can choose this 
eigenvalue to be iVj 1 . The F-term equations of the superpotential ( |2.31| ) take in this case 
the form: 



rn 
a 



M{ =0; N{ = mA5{, i,j>l; ( 2 - 38 ) 

M =—I Nc -i , 
a 

where [detiVjn is the determinant of with the first row and column discarded (i.e. the 
(11) minor of the matrix N). 

Thus, the meson matrix M takes the form 

M=rr °i (2.39) 



o o / 

where is proportional to the N c x N c identity matrix, 

M« = -I Nc . (2.40) 
a 

The (Nf - N c + 1) x (N f - N c + 1) matrix N has a block proportional to the (Nf - N c ) 
dimensional identity matrix ( |2.38|) . Overall, we find that the number of vacua in the 
magnetic gauge theory with the superpotential ( |2.4j ) agrees precisely with ( [2.14| ). 



10 



2-4- Matching the electric and magnetic descriptions 

In the previous two subsections we performed a detailed analysis of the supersymmetric 
vacua of the magnetic gauge theory (|2.4| ), and in particular reproduced ( |2.14| ) in that 
description. It is interesting to check this result in the electric description. In fact, this 
does not require any additional work. As mentioned above, the description of the electric 
gauge theory via the superpotential ( |2.6| ) is identical to the magnetic one Q2.4Q with the 
substitutions 

N c -> N f - N c , 

Am * A e , 

A-> -A , 

(2.41) 

(m, a) — > -(m e , a e ) , 
(q,q)^(Q,Q) , 
M — > AT . 

As a check, note that the transformation on the first three lines of (|2.41|) is a symmetry of 
the scale matching condition ( |2.2fl . 

To see how the vacua we found in the previous subsections arise in the electric de- 
scription ( |2.6| ), consider for example the case N c < Nf < 2N C (the regime Nf > 2N C is 
very similar). In the magnetic description we found vacua in which the meson matrix M 
( pT3| ) takes the form (up to global symmetries) 

M = diag(x Ar/_fc , y k ) (2.42) 



with x and y satisfying the relations ( |2.19| ), (|2.37[ ). The number of solutions as a function 



of k is N c - k for < k < N f - N c , and 2N C - N f for N f -N c <k<^-. 

To analyze this case in the electric variables we need to use the map (|2.41| ). Due to 



the transformation of the number of colors, the electric theory is actually in the opposite 
regime of the analysis of subsection 2.3, Nf > 2N C = 2(Nf — N c ). According to that 
analysis, the electric meson matrix takes in this case the form 



N = diag(xe f \y k e ) (2.43) 



where now for < k < N c = Nf — N c there are Nf — N c — k = N c — k solutions, and 

2 
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for Nf — N c = N c < k < ^if- there are Nf — 2N C = 2N C — Nf solutions. Comparing to 



the magnetic analysis, we see that the degeneracies are exactly the same for all k, and we 
should identify the electric and magnetic vacua for each value of k separately. 

We next show that the detailed form of the meson matrix one finds in the electric 
and magnetic descriptions is indeed the same for each k. To facilitate the comparison, we 
need to translate the results of the electric analysis, which give the auxiliary gauge singlet 
meson matrix N to those for the meson matrix M ( |1.3|) . The equation of motion of N 
arising from the superpotential (|2.6| ) gives the relation between the two: 

M 

— = m e I Nf - a e N . (2.44) 

Plugging ( 2.42|) , (|2.43| ) into ( |2.44| ) we find the following relations between the eigenvalues: 

x =A(m e — a e x e ) , 

(2.45) 

y =A(m e - a e y e ) . 



These relations can be simplified by recalling that x e and y e satisfy an analog of ( |2.19| ) 
obtained by making the replacements ( |2.41|) , 



Tfl 

x e + y e = — (2.46) 



using which we can rewrite (|2.45| ) as 

x =a e Ay e , 

(2.47) 

y =a e Ax e . 

The non-trivial check is that the polynomial equation satisfied by x, y, ( |2.37|) , and the 
corresponding equation for x e , y e , 



are compatible for all k. Plugging ( [2.471) into fl2.37j) and using the scale matching condition 



T2|) as well as fl2.7j) one finds that this is indeed the case. 

2.5. The classical limit 

In the previous subsections we found that the number of supersymmetric vacua of the 
deformed SQCD system described by (|2l| - is given by (|2TT41) . For N f > 2N C it 



agrees with the classical analysis of the magnetic theory, while in the electric description 



12 



some of the vacua exist only in the quantum theory and go off to infinity in the classical 
limit. For Nf < 2N C it is the other way around. 

To see how this happens in detail, we can take the classical limit of our general 
results. Consider, for example, the limit in which the magnetic theory becomes classical,0 
A^ Vc f — > 0. In this limit (which is equivalent via (|2.2| ) to Ae c f — > oo) the quantum 
corrections to the classical superpotential vanish (see e.g. ( |2.35|) ). To see what happens to 
the quantum vacua, we need to analyze the behavior of the eigenvalues (x, y) of subsection 
2.3 in this limit. 

Consider, for example, the case Nf < 2N C . The equations for the eigenvalues, fl2.19p , 
( |2.37|) take different forms for different k. For < k < Nf — N c one has 

= (-a) N f- N «Al N °- N f - x) N, ~ Nc ~ k . (2.49) 

This is a polynomial equation of degree N c — k. The solutions exhibit two types of behavior 
in the classical limit. There are Nf — N c — k solutions of the form 

Tft 

x= — -e, y = e, (2.50) 



with 



e N f -N c -k „ { _ a) N c -N f k A N f -SN c (2 51) 



Note that e — > in the classical limit. Thus, these solutions are small deformations of the 
classical solutions found in subsection 2.1. 

In addition, ( |2.49| ) has 2A^ C — Nf solutions in which 



x 



2N C -N f „ ^k a N f -N cA 3N c -N f 



These solutions go to x = oo in the classical limit. 



For N f -N c <k<-f- ([2T49|) is a polynomial equation of degree 2N C —Nf. Its solutions 



have the form ( j2.52| ), so they too go to infinity in the classical limit. The baryonic vacuum 
( p.38|) is a small deformation of the classical one. 

For Nf > 2N C one can perform a similar analysis and verify that all the solutions of 
( [2 .371) are small deformations of the classical ones, as expected. 



1 The analysis of the classical limit in the electric theory is very similar. 
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3. Metastable vacua 

In the previous section we discussed the supersymmetric ground states of deformed 
SQCD ( |1.4j ). In this section, we will see that this theory has many non-supersymmetric 
metastable ground states as well. As in 0, we will restrict attention to the regime Nf < 
|iV c , where the magnetic theory is free in the infrared, and can be thought of as the 
effective low energy description of the asymptotically free electric gauge theory. A similar 
analysis can be performed in the free electric phase.! 

The Kahler potential of the singlet meson superfield M in ( |2.4|) in the free magnetic 
phase is expected to take the form 

K=— 1 —M^M + --- (3.1) 

a|A e | 2 

near the origin M = 0. The positive real constant a is not easy to calculate (see for 
further discussion). It is convenient to define a superfield $ via 

M = v^A e $ (3.2) 

such that the Kahler potential for $ and the magnetic quarks is canonical near the origin 
of field space, 

K = Txq^q + Trq^q + Tr$^$ H . (3.3) 



Corrections to the Kahler potential ( p73|) are due to physics at or above the scale A m where 
the magnetic gauge theory breaks down and is replaced by the asymptotically free dual 
electric theory discussed above. The leading corrections are expected to be quartic in the 
fields and suppressed by two powers of A m , 



SK 



Tt($t$)2 _|_ . . . . (3.4) 



IA I 2 



In order to be able to ignore them we will require the expectation values of the fields q, q, 
$ to be much smaller than |A m |. 

The magnetic superpotential (|2.4|) is given by 



W mag = hq^qJ - Tr (hfiH - ^h 2 fi^ = jqtM^ + Tr Q«M 2 - mil/) . (3.5) 



2 Or, more generally, whenever one of the descriptions is weakly coupled at the energy scale 
associated with the metastable vacua. 
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By equating the two expressions in ( |3.5|) we can relate the parameters as follows: 

h = 5 A* 2 = m ^ ■> = a ^- 2 ■ (3-6) 



The classical supersymmetric vacua (|2.10|) , ( |2.11|) take in these variables the form 



-=(" 2 o 4 o)- ^ 

To construct the metastable states, it is convenient to further split the (Nf —k)x (Nf — k) 
block at the lower right corner of (|3.7| ), (|3.8| ) into blocks of size n and Nf — k — n as follows: 

/0 \ 

h$ = h$ n \ (3.9) 

V° fa*,-*-*/ 

and 

/V4 0\ 
qq = j (pip . (3.10) 
\ 0/ 

and </5 are n x (A^j — N c — k) dimensional matrices. They correspond to n flavors of 
fundamentals of the gauge group SU(Nf — N c — k) which is unbroken by the non-zero 
expectation value of q, q in ( |3.10j ). $ n and (pup are n x n matrices. The supersymmetric 

2 

ground state ( p.7|) , (|3.8|) corresponds to /i$ n = ^-I n -, <p = <p> = 0. 

As we will see next, there are metastable vacua near the origin as well. We will restrict 
to the regime 

A m >/i>/^ (3.11) 

where the first inequality is as in 0, and the second implies that the term proportional 
to fXff, is a small perturbation of the superpotential considered in 0. Since in general the 
expectation value of $ ( |3.9| ) can be large in the regime ( |3.11fl , it is convenient to discuss 
separately the cases n = Nf — k and n < Nf — k. 
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3.1. Metastable vacua with n = Nf — k 

In this case, the dynamics of the fields $ n , (p, <p near the origin of field space is 
obtained from the underlying SU(Nf — N c ) gauge theory by giving an expectation value! 
fitok flavors of magnetic quarks q, q, ( |3.1U| ) . The low energy theory is an SU (Nf — N c — k) 
gauge theory with Nf — k = n light flavors, which is infrared free in the regime Nf < %N C . 
Its scale, A;, is related to that of the underlying theory, A m , via the scale matching relation 



A useful way of writing ( |3.12| ) is 



\ 2k /. ^ 3(N f -N c -k)-n 



Am J V A 



(3.13) 



In the regime ( |3.11| ) , the left hand side of (|3.13|) is very small. Since the power on the right 



hand side is negative, we conclude that A m ^> A;. Using the fact that ( |3.12| ) also implies 
that 



in 

AJ \A~ 
we conclude that the hierarchy of scales is 



x2fc , K \ 3(Nf-N c )-Nf 



(3.14) 



H < Ai < A m . (3.15) 

Thus, we see that as long as fx A m , the gauge dynamics is weakly coupled for all 
energies well below A m . At energies above fi, the full magnetic gauge group is restored. At 
an energy of order \i the theory crosses over to the one with gauge group SU (Nf — Nc — k), 
which is also infrared free and due to ( |3.15| ) is weakly coupled. Therefore, we can mostly 



neglect the gauge dynamics in what follows. 

The potential for $ n , ip, (p near the origin of field space has two relevant contributions. 
One is the tree level potential that follows from ( [373]) , (|3.5| ). The second is the one loop 



potential, which to leading order in h, ^ is identical to that computed in J7[ , and is given 
by 

y 1 _ loop = 6|/iV| 2 Tr$J l $ n (3.16) 



3 



In the regime ( |3.11 ), the expectation value of q is such that the corrections to the Kahler 



potential discussed around eq. (3.4) can be neglected, as in 0. 
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with b a numerical constant, 



In 4 - 1 

8tt 2 



(JV> - N c ) . 



The full one loop potential for $ n , ip, tp has the form 

! : |$ n ^| 2 + |<^$ n | 2 + \<p<p - n 2 I n + h^ n \ 2 + b\h^\ 2 Tr^ n $ r 



h 



(3.17) 



(3.18) 



We are looking for a minimum of the potential at tp 
we find that if there is one, it is located at 



tp = 0. Differentiating (|3.18|) w.r.t. 



2 J ' n ' 



and its vacuum energy is given to leading order in h, — by 



2 1 2 



(3.19) 



(3.20) 



Expanding around this solution one finds that the mass matrix of f, if has eigenvalues 



in 



1 



± 



(l/^l 2 + % 



2^2 



[|^| 2 ± b\h\ 2 {\^\ 2 + b\^\ 2 )} ~ - (|^| 2 ± |6V| 2 ) . (3.21 



Equation ( |3.21 ) implies that for sufficiently small some of the fluctuations of ip, (p are 
tachyonic. To avoid tachyons one must have 



> 



l-6|7i| : 



(3.22) 



The condition ( |3.22| ) is compatible with \h\, ^ being arbitrarily small. If it is satisfied 



the quarks tp, <p are massive, and the vacuum (|3.19| ) is locally stable. 

For k < Nf — N c — 1, this vacuum contains an unbroken SU(Nf — N c — k) gauge theory, 
which is weakly coupled at energies above the scale of the masses of (p, <p (|3.21|) . For energies 
well below these masses this gauge theory confines and has, as before, Nf — N c — k vacua. 
For k = Nf — N c — 1 there are no unbroken gauge fields and the theory is weakly coupled 
at long distances. 

For k = Nf — N c (the largest value it can take), the fields cp, <p in ( |3.10| ) do not exist. 
Hence, in this case the constraint ( |3.22|) is absent and the resulting metastable vacuum 
is present for arbitrarily small In the limit \x§ — > it corresponds to the metastable 
state discussed in 0. Even when ( |3.22[ ) is satisfied, the vacuum with k = Nf — N c is more 
long-lived than those with k < Nf — N c since the mode of instability towards condensation 
of <p, <p is absent in this case. 
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3.2. Metastable vacua with n < Nf — k 

In order to trust an analysis based on the Kahler potential (|3.3|) we have to demand 
that all components of $ ( |3.9|) are much smaller than A m . In the present case this implies 

u 2 

— < hA m . (3.23) 



A useful way of thinking about (|3.23|) is as the requirement 



« , (3.24) 

Am /" 

which is compatible with all three couplings and — being small. Furthermore, this 

requirement is very natural in the brane realization of the theory, |3l] 



Like in the previous subsection, we would like the gauge dynamics to be weak near 
the origin of the field space of (p, <p, $ n , ( |3.9| ), ( |3.10| ). The dynamics of these fields is 

2 

obtained from the magnetic SU(Nf — N c ) gauge theory by giving masses ^ to Nf — k — n 
flavors (see (|3.9| )) and vacuum expectation values \x to k flavors (see ( p.lOp ). In the regime 
( p.ll|) , the masses are much larger than the expectation values. Thus we can analyze the 
reduction from the SU(Nf — N c ) magnetic gauge theory to the low energy SU(Nf — A^ c — k) 
one in two steps, by first incorporating the masses, and then the expectation values. We 
would like the theory to remain weakly coupled throughout this process. 

In the first step, we go from an SU(Nf — N c ) SYM theory with Nf flavors and scale 
A m to one with the same gauge group, n + k flavors and scale Ai. The scale matching 
relation between the two theories is 

(2 \ Nf—k—n 
= A 3(^-^)-(fc+n) _ (325) 

The fact that the magnetic gauge theory is not asymptotically free is reflected in the power 
of A m on the left hand side of (|3.25| ) being negative. The theory with n + k flavors can 
be either asymptotically free or not, which is reflected in the fact that the power of Ai on 
the right hand side can be positive or negative. 
Consider first the case@ 

k + n > 3(N f - N c ) (3.26) 



4 In the free magnetic phase Nf > 3(Nf — N c ) there are solutions to this constraint. For large 
generic Nf, N c , the number of such solutions is large as well. 
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in which it is not asymptotically free. Rewriting ( 3.25 ) in the form 

Nf-k-n 

^ N 3(N f -N c )-k-n 
Am. 




(3.27) 



and using the fact that the left hand side is very small due to (|3.23| ), we conclude that 

Ai > A m . (3.28) 
Hence the SU(Nf — N c ) gauge dynamics is weakly coupled both for energies larger than 

2 

— , for which the number of flavors is Nf and for lower energies, where it is n + k. 

In the second step of the reduction we give an expectation value fx to k flavors of 
quarks and go down to an SU(Nf — N c — k) gauge theory with n light flavors and scale 
A;. The scale matching relation is (compare to ( |3.12| )) 

^(N f -N c )-(k+n) _ ^2k^3(N f -N c -k)-n (3 29) 

In the regime (|3.26|) the power of A; in (|3.29|) is negative. Thus, the low energy SU(Nf — 
N c — k) gauge theory with n flavors is weakly coupled below the scale A;. Using the fact 
that (i <C Ai, one can see from Q3.29|) that A; is in the range 

H < Ai < Ax . (3.30) 

We see that the theory is weakly coupled at all scales. At energies much above \i it is 
governed by the intermediate theory with scale Ai, which as we saw before is weakly 
coupled. At the energy \x it crosses over to the low energy theory with the scale A;, 
which due to ( |3.3U[ ) is weakly coupled as well. Consequently, the gauge dynamics can 
be neglected, as in the previous subsection, and one can proceed as there, with similar 
conclusions. 

To summarize, we find that deformed SQCD in the regime ( |3.11|) , ( |3.24j ) in coupling 
space has metastable vacua of the form (|3.9|) , ( |3.10| ), (|3.19|) , with k and n satisfying the 
constraint (|3.26| ) . The states with k < Nf — N c further require (|3.22|) for their existence. 
Those with k = Nf — N c exist throughout this regime. 

The above analysis can be repeated for the case where the constraint ( |3.26| ) is not 
satisfied. One finds that there are regions in parameter space in which the gauge dynamics 
is weakly coupled and metastable states exist. We will not describe the details here. 

Finally, it is interesting to ask what happens to the states with k = Nf — N c , n < N c as 
Htjj — » 0. As decreases, at some point we leave the weak coupling regime ( |3.23|) and the 
analysis presented above becomes unreliable. The limit — >• involves strongly coupled 
dynamics. The brane construction of [31]] seems to suggest that these vacua survive in the 
limit, and thus in massive SQCD there are additional metastable states to those studied 
in R]. More work is required to resolve this issue conclusively. 
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4. Discussion 



In this paper we discussed supersymmetric QCD in the presence of the superpotential 
(|2.3|) . We described the supersymmetric vacua of the theory, and generalized the discussion 
of metastable supersymmetry breaking states in [[7] to this case. We saw that by tuning 
the parameters of the theory one can make the non-supersymmetric states arbitrarily long- 
lived. In other regions of parameter space, which can be studied using weakly coupled field 
theory, many of these states become unstable and disappear. 

The mechanism for the emergence of metastable states in the theory with a 7^ 
Q2.3|) is slightly different than that in the theory with a = studied in 0]. There, the 
classical theory spontaneously broke supersymmetry and had a pseudo-moduli space of 
non-supersymmetric states. One loop effects lifted the pseudo-moduli space and replaced 
it by an isolated metastable state. In our case, classically there are supersymmetric ground 
states and no non-supersymmetric metastable states. The latter are due to a competition 
between the classical and one loop contributions to the potential. The fact that one can 
ignore higher loop contributions to the potential in studying these states is analogous to 
what happens in the e expansion in quantum field theory. 

The theory with quartic superpotential in (|2.3| ) requires a UV completion. One possi- 
bility is to introduce the singlet meson N, and rewrite it as ( |2.6|) , which is renormalizable. 
Alternatively, one can proceed as follows. Consider first SQCD with vanishing superpo- 
tential and Nf < 3N C . Due to asymptotic freedom, it approaches a free fixed point in the 
UV, while at long distances it is governed by a non-trivial fixed point, at which the scaling 
dimension of the meson operator (|1.3|) is smaller than its free field value. For Nf < 2N C , 
the IR dimension is sufficiently small that the quartic perturbation in (2.3) becomes a rel- 
evant perturbation of this fixed point. Thus, in this regime the theory with superpotential 
( |2.3| ) is well defined in the UV. Our analysis of metastable states was performed (following 
0) in the region Nf < |iV c , in which the magnetic theory is infrared free and the electric 
theory is UV complete. 

A natural generalization of the problem studied here is to higher order superpotentials 
of the form ( |1.2| ). The supersymmetric ground states can be analyzed as in section 2. For 
example, in the mesonic branch, in which the matrix M is non-degenerate, the F-term 
equation ( |2.36|) takes in the general case the form 



1 / del M ^ 
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Since the left hand side is a polynomial of degree no, the matrix M has at most no distinct 
eigenvalues, (x\, X2, ■ ■ • , x no ). In a general vacuum, it takes the form 

M = di a g(x[\x l 2 \---,xn °) (4.2) 

with lj = 0, 1, ■ ■ • , Nf and ^ ■ lj = Nf. It seems clear that for large N the number of 
vacua grows like N n °, where N denotes collectively Nf, N c and any linear combinations 
of the two. The growth of metastable states with N is even faster since we can follow the 
construction of section 3 in each block of M ( |4.2|) separately. 

From the perspective of effective field theory it is natural to consider superpotentials 
with no of order Nf, N c . The number of stable and metastable vacua grows in this case 
at least as fast as N\. To ensure that these vacua are all long-lived one needs to fine 
tune the couplings m n in ( fl.l|) accordingly. Conversely, if the couplings are generic, one 
expects a distribution of lifetimes, so that at least some of the non-supersymmetric vacua 
are long-lived. 

An interesting]! generalization of the analysis of this paper is to a system in which the 
SU(Nf)<img symmetry of ( |1.2|) is gauged. This leads to adjoint SQCD, which was studied 
using a generalization of Seiberg duality in [p5|-p7|. One can use these results as well as 



those of this paper to analyze the metastable states in this system. 

As mentioned in the introduction, SQCD with the superpotential (|2.3| ) has a simple 
embedding in string theory, as a low energy theory on a system of intersecting D-branes 
and A^S'5-branes. In a companion paper |31| we describe this embedding and show that 



much of the structure found in the gauge theory description in this paper is nicely realized 
in the brane system. 
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